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Within the framework of stochastic inflationary cosmology we derive steady-state distributions 
P C {V) of domains in comoving coordinates, under the assumption of slow-rolling and for two specific 
choices of the coarse-grained inflaton potential U(<3>). We model the process as a Starobinsky-like 
equation in U-space plus a time-independent source term P W (V) which carries (phenomenologically) 
quantum-mechanical information drawn from either of two known solutions of the Wheeler-De Witt 
equation: Hartle-Hawking's and Vilenkin's wave functions. The presence of the source term leads 
to the existence of nontrivial steady-state distributions P™(V). The relative efficiencies of both 
0^ ■ mechanisms at different scales are compared for the proposed potentials. 

On 
> 

Since the differential microwave radiometer (DMR) mounted on the COBE satellite first detected temperature 
anisotropies in the cosmic background radiation, we have the possibility to directly probe the initial density pertur- 
bation. The fact that the resulting energy density fluctuations (-^ w .5x 10~ 5 ) fit the scaling spectrum predicted by 
the inflation model, suggests that they had indeed their origin in the quantum fluctuations of the "inflaton" scalar 
field during the inflationary era. Although this problem is in principle of a quantum-mechanical nature, the fact that 
l ' under certain conditions (which are made precise in Ref. 0) the inflaton field can be considered as classical largely 
. simplifies the approach, by allowing a Langevin-Ytke stochastic treatment. Stochastic Inflation branched off and grew 
up from within the context of chaotic inflation 0. It was Starobinsky || the first one to derive, from a Langevin- 
likc equation for the dynamics of the coarse-grained inflaton field <& and in the absence of back-reaction from the 
metric, a Fokker-Planck equation for the transition probabilities P c ($, i') in comoving coordinates. The latter 
provide us with statistical information about the relative number of spatial domains that evolve in a time interval 
t — t' from having a typical value (assumed constant throughout the domain) towards a new configuration with a 
| typical value $. This approach is purely phenomenological and keeps very little information about its quantum origin, 
— «. ■ since the coarse-graining procedure erases quantum correlations. Hence, the issue of integrating quantum- mechanical 
^ i information to the inflationary evolution is still a pending assignment. 

I In a recent paper Q we have offered a phenomenologically drawn example that takes into account — if not properly 

back-reaction — at least the possibility that Hubble's constant H depend on Our approach disregarded instead the 
boundary conditions at Planck's and exit-of-inflation $ regimes and, in particular, domain injection from the purely 
quantum-mechanical phase. As we know, the Wheeler-Dc Witt (WDW) equation is a possible way to describe the 
k>( | fully quantum-mechanical behaviour of the Universe. However, this equation has infinitely many solutions: hence, 
there have been several attempts to solve it by starting from different initial conditions, each providing a different 
mechanism whereby domains with a typical size a = H^ 1 can enter the inflationary stage. Two very popular ones 
among them are quantum creation from "nothing" Q and the no boundary proposal 

It is our aim in this work to investigate the inflationary stage in the light of the information provided by solutions 
to the WDW equation. Our (still phenomenological) approach will be to enlarge the scope of Starobinsky's equation 
by adding a time-independent source term. We shall take this (extrinsic) domain-injection probability per unit time 
to be proportional to the squared modulus of the aforementioned solutions. Moreover, we shall assume that during 
the inflationary phase the slow-rolling condition is satisfied. Under these assumptions, we obtain phenomenologically 
meaningful stationary domain distributions in U-space for two physically appealing forms of the inflaton potential 

For the role that quantum cosmology is intended to play in this work (i.e. to act as an effective source of homogeneous 
horizon-sized domains) it suffices to consider the minisuper space version of the WDW equation Q: 
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This is the cornerstone for a quantum-mechanical description of the pre-inflationary era, and is derived under the 
assumption that the metric is of the Friedmann-Robertson- Walker (FRW) type 



ds 2 = a 1 [N 2 {t)dt 2 ~ a 2 {t)dttl] 



(2) 



and that the only quantum variables in this stage are the FRW scale factor a and the inflaton field $. In eq.(g) 
N(t) is the lapse function, dil 2 the metric on a unit three-sphere and a 2 = 2G/3n a normalizing factor chosen for 
convenience. Equation (|l]) arises as a constraint when varying the Hamiltonian with respect to N. 

The asymptotic solution of eq. (jl|) with the boundary condition \P(a — * oo) — > was obtained by Vilenkin j^j, in 
analogy with the decay of a metastable state. Taking a = ii" _1 (<i>), its squared modulus 



P V (V) = exp (-3M p 4 /8y) 



(3) 



describes the nucleation of universes (domains in our phenomenological model) with a given value of the inflaton 
potential V in the inflationary phase through quantum tunneling, and the corresponding process is known as "quantum 
creation from nothing" . On the other hand, the celebrated Hartle-Hawking (HH) solution M leads to 



P h (V) = exp {3M*/8V). 



(4) 



We find that eq.(^) and eq.(Q) differ in the sign of the exponential. In the following — and unless otherwise specified — 
we shall denote the domain-injection probability arising from any solution to the WDW equation as P W (V). Obviously, 
given a specific functional form V(&) one immediately obtains P tu ($). 

Now we want to incorporate this information into the equation for P C {V, t) (the probability distribution of domains 
in comoving coordinates during the inflationary stage). After choosing units such that M p = 1, changing variables in 
the Fokker-Planck equation for P c (<E>,i) and recalling that = V -§y (the prime denotes ^), one obtains 



V^-^V'— (viPp* 
dV 
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+ C W P W (V), 



where the time-independent domain distribution P W (V) acts as a source term that accounts for domain injection into 
the inflationary phase, and C w is a (unknown) constant. /3 = 1/2 for the Stratonovich prescription and (3 = for 
Ito's one. 

In this work we shall restrict ourselves to the stationary equation which, by defining 



g w (v) = v-^-Vv'^ (yWp?) 
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can be rewritten as 
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with Ci — %\/QttC w , and G^ l (V) is the function of V solution of the last equation. Once we solve eq. 
obtain P?(V) from eq.(g), setting (3 = 1 as in Ref. @: 



P^(V) = K V-^ 2 e Xp [ w 
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Here B is an integration constant and K Q = V~ x — 1 (with V e — V(<j> e )) comes from normalizing at the exit-of-inflation 
value <& e of the inflaton field. 

In the remaining of this letter — and with the aim to illustrate the procedure — we derive the steady-state domain 
distributions P%(V) and P^{V) arising from both domain-injection mechanisms, for two specific choices of V(3>). 



1. Let us first consider the case of a quadratic potential V($>) = 1 ^-^> 2 , so that V' = (2m 2 V) 



1/2. 
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(a) If the domains are produced by quantum tunneling, P W (V) is Vilenkin's distribution eq.(|^). For simplicity 
we choose Cf — — 1. This gives 
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where Erf is the error function. The general solution of eq.(Q) is 
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P*(V) = K„ exp ( — 



(8) 



where Bj, is a constant of integration and 
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From the requirement that P^(V) be positive in the inflationary regime it follows that > | min[/i(V)]|. 
We need another condition in order to fully specify the integration constants. Considering our ignorance, 
we will require both solutions to agree at the Planck scale i.e. that P% (V p ) — P V {V P ), where V p = V($> p ). 
Then 
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(b) In the case of a HH source term [|| and choosing C\ = 1, it results 
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The steady-state domain distribution in comoving coordinates is then 



P?(V)=K expl^-) bP- 9i (V) 



where 
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and is an integration constant required now to be B^' > \ min\gi(V)]\. Imposing as before P^(V P ) — 



Pv(V p ), 



B 
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The functions fi(V) and gi(V) — describing the domain flow from the quantum sector towards the inflationary 
regime — must be evaluated numerically in this example. 

2. The following is an exactly solvable example: 

(a) Let us consider in equation (||) the choice P V (V)/V = V~ 2 , and choose also Cf = —1. The resulting 
scalar potential is 
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The condition V($) > imposes the restriction $ p — $ > 8/3. Since, on the other hand G s v {V) 
it is possible to find P^ (V) from equation (Q): 
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Here By' is again a constant ensuring that the distribution P£{V) be positive in the whole inflationary 
domain. By choosing P£ (V p ) — P V (V P ) and taking V p = 1, it results to be 
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(b) When the source term is HH-like, choosing P h (V)/V = exp (^)— so that = § exp (^)— the 

solution for P^ (with C{ 1 = 1) results: 



P e »(V) = A-.V- 3/2 exp( — 
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By choosing P c "(F p ) = P V (V P ), 
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Figure 1 compares the ratio Z(V) = P^(y)/P c ' I (V r ) as a function of V for the two cases, between the values V e » 
and VJ, = 1. Although both curves are rather similar in form when looked in their own scale, Z(V) for the logarithmic 
potential looks very flat when compared with the quadratic one. Moreover — except for V < 0.1 — the quadratic 
potential predicts larger values of Z(V). In fact, whereas for the quadratic potential and V > 0.1 it is Z(V) > 1 (with 
a maximum larger than 3 at V ~ 0.3) — thus implying a higher efficiency of the nucleation mechanism as compared to 
the HH one — for the logarithmic one it is always Z(V) < 1 (thus implying a higher efficiency of the HH mechanism) 
and this relative efficiency passes through a (less dramatic) minimum at V ~ 0.2. 

Summing up, in this work we have studied the effect of the domain-injection mechanism on the (stationary) domain 
distribution at the end of the inflationary era, for two phenomenologically appealing proposals of inflaton potential 
forms: quadratic and logarithmic, respectively. The domain-injection mechanisms are provided in our approach by 
known solutions to the WDW equation. Although the ratio Z(V) remains of order one and looks similar for both 
potential forms, it is concluded that the relative efficiency of each mechanism depends on the form of the inflaton 
potential. 

This letter is a first step towards studying the possible appearance of stability in Stochastic Inflation due to quantum 
sources of inflationary domains, deserving the matter further research and debate. We are aware of the lack of natural- 
ness of our approach, that combines such different descriptions as quantum cosmology and (gravitationally classical) 
stochastic inflation. The transition from quantum to classical behaviour in quantum cosmology, and particularly 
the conditions under which the wave function of the Universe becomes semiclassical, are far from being understood. 
Thus, considering the possibility that the exit of the quantum phase is not a well defined event, a stochastic interplay 
between both classical and quantum descriptions seems plausible. 
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FIG. 1. The ratio Z(V) = P"(V)/P^(V) between the steady-state distributions arising from Vilenkin's and Hartle-Hawking's 
source terms, for the case of a quadratic potential (upper curve) and for a logarithmic one (lower curve) . 
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